STATIONARY DISTRIBUTIONS FOR JUMP PROCESSES WITH 

INERT DRIFT 



K. BURDZY T. KULCZYCKI AND R.L. SCHILLING 

Abstract. We analyze jump processes Z with "inert drift" determined by a "mem- 
ory" process S. The state space of (Z, S) is the Cartesian product of the unit circle 
and the real line. We prove that the stationary distribution of (Z, S) is the product 
of the uniform probability measure and a Gaussian distribution. 



1. Introduction 

We are going to find stationary distributions for jump processes with inert drift. We 
will first review various sources of inspiration for this project, related models and results. 
Then we will discuss some technical aspects of the paper that may have independent 
interest. 

This paper is concerned with the following system of stochastic differential equations 
( the precise statement is in the next section), 

(1.1) dY t = dX t + W'(Y t )S t dt, 

(1.2) dS t = W"(Y t )dt, 

where X is a stable Levy process and W is a C 5 function. This equation is similar to 
equation [TJ (4.1)], driven by Brownian motion, but in (11.1 ft the term ~ (AVV)(X t ) dt 
from the first line of [1, (4.1)] is missing. An explanation for this can be found in 
heuristic calculations in [SJ Example 3.7]. The paper [S] deals with Markov processes 
with finite state spaces and (continuous-space) inert drifts. This class of processes is 
relatively easy to analyze from the technical point of view. It can be used to generate 
conjectures, for example, [8j Example 3.7] contains a conjecture about the process 
defined by (O^ - flO]) . 

The main result of this paper, i.e. Theorem I2.12[ is concerned with the stationary 
distribution of a transformation of (Y, S). In order to obtain non-trivial results, we 
"wrap" Y on the unit circle, so that the state space for the transformed process Y is 
compact. In other words, we consider (Z t ,S t ) = (e lY *,S t ). The stationary distribution 
for (Z t , S t ) is the product of the uniform distribution on the circle and the normal 
distribution. 

The product form of the stationary distribution for a two-component Markov process 
is obvious if the two components are independent Markov processes. The product form 
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is far from obvious if the components are not independent but it does appear in a 
number of contexts, from queuing theory to mathematical physics. The paper [8] was 
an attempt to understand this phenomenon for a class of models. 

One expects to encounter a Gaussian distribution as (a part of) the stationary dis- 
tribution in some well understood situations. First, Gaussian distributions arise in the 
context of the Central Limit Theorem (CLT) and continuous limits of CLT-based mod- 
els. Another class of examples of processes with Gaussian stationary measures comes 
from mathematical physics. The Gibbs measure is given by C\ exp(— c<i Ylij( x i ~ x j) 2 ) 
in some models, such as the Gaussian free field, see [TTj . In such models, the Gaussian 
nature of the stationary measure arises because the strength of the potential between 
two elements of the system is proportional to their "distance" (as in Hooke's law for 
springs) and, therefore, the potential energy is proportional to the square of the distance 
between two elements. Our model is different in that the square in the exponential func- 
tion represents the "kinetic energy" (square of the drift magnitude) and not potential 
energy of a force. The unexpected appearance of the Gaussian distribution in some 
stationary measures was noticed in [7] before it was explored more deeply in [Sid]. 

The present article has a companion [5] in which we analyze a related jump process 
with "memory" . In that model, the memory process affects the rate of jumps but it 
does not add a drift to the jump process. The stationary distribution for that model is 
also the product of uniform probability measure and a Gaussian distribution. 

An ongoing research project of one of the authors is concerned with Markov processes 
with inert drift when the noise (represented by X in fll.il) ) goes to 0. In other words, 
one can regard the process (Y, S) as a trajectory of a dynamical system perturbed by a 
small noise. No matter how small the noise is, the second component of the stationary 
measure will always be Gaussian. Although we do not study small noise asymptotics 
in this paper, it is clear from our results that the Gaussian character of the stationary 
distribution for the perturbed dynamical system does not depend on the Gaussian 
character of the noise — it holds for the stable noise. 

Models of Markov processes with inert drift can represent the motion of an inert 
particle in a potential, with small noise perturbing the motion. Although such models 
are related to the Langevin equation (see [13]), they are different. There are several 
recent papers devoted to similar models, see, e.g., [21 [3j HJ |5] . 

We turn to the technical aspects of the paper. The biggest effort is directed at 
determining a core of the generator of the process. This is done by showing that the 
semigroup T t of the process (Y t , S t ) preserves C%, see Theorem 13.11 The main idea is 
based on an estimate of the smoothness of the stochastic flow of solutions to f ll.ip - 
(11. 2p . This result, proved in greater generality than that needed for our main results, is 
presented in Section [31 see Proposition [373] This proposition actually makes an assertion 
on the pathwise smoothness of the flow. It seems that Theorem 13.11 and Proposition 
13.31 are of independent interest. 

1.1. Notation. Since the paper uses a large amount of notation, we collect most fre- 
quently used symbols in the table below, for easy reference. 
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a V b, a A b 


max(a, &), min(a, b); 


a+, a_ 


max(a, 0), — min(a, 0); 




m 

E Xjl where x = (xi, . . . , x m ) G R m ; 


efe 


the fc-th unit base vector in the usual orthonormal basis for R"; 




/ 1 -t- rv\ 1 a ~^ 

" r ( 2 )^r(i-«- 




, a = (ai , . . . , cm) e Wn; 

dx" 1 ■ ■ ■ dxj d 


c fe 


fc-times continuously differentiable functions; 


/~ifc /~ik /-~ik 


functions in C k which, together with all their derivatives up to 
order k, are "bounded", are "compactly supported", and "vanish 
at infinity", respectively; 


ll/lloo.B 


sup \f(x)\ for / : R™ -> R; 


P 0) /l|oo, S 


E P'VIU*; 

\a\=3 


II/II0),b, H/Hci) 


E Bup|£> a /(*)|,resp- E H^/IU 

| Q |<i xeB |a|<j 




77. 71 

sup D a Vfc(x)|, resp., ^ ^ D Q T4 oo for any function 

\a\—j k—1 | a —j k—1 

V : R" -S- R"; 




j J 

Ell^ V H-,B,resp.,EP W ^lloo; 

i=0 i=0 


s 


{zeC : |z| = 1} unit circle in C. 



Constants c without sub- or superscript are generic and may change their value from 
line to line. 



2. A JUMP PROCESS WITH A SMOOTH DRIFT 

Let S = {z G C : \z\ = 1} be the unit circle in C. Consider a C 5 function V : S ->■ R 
which is not identically constant and put W(x) = V(e tx ), x G R. Let X t be a symmetric 
a-stable Levy process on R which has the jump density A a \x — y\~ l ~ a , a G (0, 2). Let 
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(Y, S) be a Markov process with the state space R 2 satisfying the following SDE, 

dY t = dX t + W'(Y t )S t dt, 
dS t = W"(Y t )dt. 



(2.1) 



Lemma 2.1. The SDE (12. ip has a unique strong solution which is a strong Markov 
process with cadlag paths. 



(2.2) 



Proof. For every n G IN define the function f n : R — > R by f n (s) := (—n) V s An. We 
consider for fixed n G IN the following SDE 

'dY t {n) =dX t + W'(Y t {n) )f n (S { t n) )dt, 

dSl n) = W"(Y t {n) )dt. 

Note that R 2 E3 (y, s) i— > W'(y)f n (s) is a Lipschitz function. By [TJJ Theorem V.7] and 
[TU Theorems V.31, V.32] the SDE (12.21) has a unique strong solution which has the 
strong Markov property and cadlag paths for every fixed n G IN. 

Now fix t < oo and a starting point R 2 3 (y,s) = (Yq , Sq^). Note that for any 
t < to we have 

r-t 



(n) 



< \s\+t \\W"\ 



Pick n > \s\ + to||W^"||oo; n G IN. For such n and any t < t , the process (Y t , S t ) := 
(Yt, St) is a solution to (12. ip with starting point (y, s). This shows that for any fixed 
starting point (y, s) = (Y , S ) and fixed t < oo the SDE (12. ip has a unique strong 
solution up to time to- The solution is strong Markov and has cadlag paths. Since 
to < oo and the starting point (y, s) are arbitrary, the lemma follows. □ 

We will now introduce some notation. Let IN be the positive integers and INo = 
IN U {0}. For any / : S ->• R we set 

/(■'•): fU u ). XGR. 

We say that / : S — > R is differentiable at z = e lx , x G R, if and only if / is differentiable 
at x and we put 

f\z) := (f)'(x), where z = e ix , x G R. 

Analogously, we say that / : S — > R is n times differentiable at z = e lx , x G R, if and 
only if / is n times differentiable at x and we write 

= (/) (n) (x), where z = e lx , x G R. 

In a similar way we define for / : S x R — y R 

(2.3) f{v,s) = f{e i »,s), y,seR. 

We say that D a f(z, s), z = e w , y, s G R, a G IN 2 ,, exists if and only if D a f(y, s) exists 
and we set 

D a f(z,s) = D a f(y,s), where z = e iy } y, s G R. 
When writing C 2 (S), C 2 (S x R), etc., we are referring to the derivatives defined above. 
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Let 

(2.4) Z t = e iYt . 

Then (Z, S) is "a symmetric a-stable process with inert drift wrapped on the unit 
circle" . In general, a function of a (strong) Markov process is not any longer a Markov 
process. We will show that the "wrapped" process (Z t , S t ) = (e tYt ,S t ) is a strong 
Markov process because the function W(x) = V(e lx ) is periodic. 

Lemma 2.2. Let (Y t , S t ) be the solution of the SDE (EH]) . Then 

F (y+2n,s)(y t g ^ _|_ 2lT , St €l B) = P fe ' S) (F t G A, S t G B) 

holds for all (y, s) G R 2 and all Borel sets A,BcE. 

Proof. Denote by (Y^, Sf) the unique solution of the SDE (12. ip with initial value 
(Yq, Sq) = (y, s). We assume without loss of generality that X = 0. By definition, the 
process (F/ +2,r , S^) solves 

y + 2ix + X t + [ W'{Y r )S r dr, 
Jo 

W"(Y r ) dr. 



Since the function W is periodic with period 2ir, we know that W'(Y r ) = W'(Y r — 27r) 
and W"(Y r ) = W"(Y r - 2tt). Therefore, (Y t y+2n , Si) solves the system 

= y + 27T + X t + [ W'(Y r -2ir)S r dr } 

= s+ [ W"{Y r - 2tt) dr. 
Jo 

By subtracting 2ir from both sides of the first equation we get 

Y t - 2tt = y + X t + J W'{Y r - 2it)S r dr, 

S t = s + W"(Y r - 2tt) dr. 

Since the solutions are unique, this shows that (l^ y+27r , S t ) = (Y? + 2tt, S t ) from which 
the claim follows. □ 

We can now use a rather general result on transformations of the state space due to 
Dynkin [9l 10.25, Theorem 10.13], see also Glover [11] and Sharpe (161 Section 13]. 

Corollary 2.3. Let 7 : R 2 — > S x R ; 7(2/, s) := (e iy , s) and (Y t , S t ) be the unique, 
cddldg strong Markov solution of the SDE (12. ip . Then (Z t ,S t ) = (e lYt ,S t ) is also a 
strong Markov process. Let Pt((y,s),A x B) denote the transition function of (Y,S) 
and P®((y, s), A x B) the transition function of (Z, S). Then for y, s G R and Borel 
sets A,BcR, 

P, s ( 7 (y, s), A x B) = P t {(y, slj-^A x B)) 
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Proof. All we have to do is to verify Dynkin's condition [9, 10. 25. A] saying that 

P t ((y, s), T\A x B)) = P t {{y\ s'),'y~ 1 (A x B)) 

holds for all Borel sets A C S, B C R and all points (y , s) , (y' , s') G R 2 such that 
7(y, s) = j(y', s'). Clearly s = s' and y — y' = 2jir for some j G Z. Denote /(y) = e iy 
. Applying Lemma I2.2I repeatedly we find 

{{Y t , S t ) G 7 _1 (A x B)) = (Y t G r\A), S t G B) 

= F iy+2nj ' s) (Y t G + 27Tj, $ G B) 

= P^+ 2 ^» (K t G 5 t G B) 

= p^+ 2 ^» ((y t , s t ) g 7 -1 (^4 xB)). □ 

We are going to calculate the generators of the processes X t , (Y t , S t ) and (Z t , S t ). 

By 9 X let us denote the generator of the semigroup, defined on the Banach space 
(Cb(IR), || ■ Hoc), of the process X t . By D(S X ) we denote the domain of S x ■ It is well 
known that C 2 (R) C D(5 X ) and for / G C 2 (R) we have 5 X f = -(-A) a / 2 /, where 

-(-A) a / 2 f(x) = A a lim / f l V) ~!!l] dy, x G R. 

If / G C 2 (R) is periodic with period 27r then we have 

-(-^ a/2 m = A a inn / ^ (y) ~fff 

(2 ' 5) „ /" f(y)-f(x) , 

^ K>\v-*\ \x-y + 2mr\i+" 



In the sequel we will need the following auxiliary notation 

Definition 2.4. 

C^R 2 ) := {/ : R 2 ->■ R : 3N > supp(/) C R x [-N, N], 

f is bounded and uniformly continuous on R 2 }, 
C 2 (R 2 ) :=C,(R 2 )nC 2 (R 2 ). 

Let us define the transition semigroup {T t }t>o of the process (Y t , S t ) by 

(2.6) T t f(y,s) = E^f(Y t ,St), y,seR, 

for functions / G C fe (R 2 ). Let 5 {Y ' S) be the generator of {Tj t > and let D(9 {Y ' S) ) be 
the domain of S <T ' S ' ) - 

Lemma 2.5. We have C 2 (R 2 ) C ^(9^) and for f G C 2 (R 2 ), 

(2.7) S (y ' 5) /(l/, s) = -(-^y) a/2 f(y, s) + W(v)«/y(v, s) + W"(y)f s (y, s), y, s G R. 
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Proof. Let / G C*(R 2 ). Throughout the proof we will assume that supp(/) C R x 
(— M ,M ) for some M > 0. Note that for any starting point (Y , So) = (y,s) G 
R x [-M , M ] and all < t < 1, 

rt 



\St\ 



S + [ W"{Y r )dr 
Jo 



<M + ||W"| 



Put 



M 1 = M + H^lco. 

Note that if (y,s) £ R x [-Mi, Mi] and (Y , S ) = (y,s) then for any < t < 1 we 



have 



\St\ 



S + [ W"{Y r )dr 
Jo 



> M 1 - \\W"\ 



M , 



and, therefore, f(Y t ,S t ) = 0. It follows that for any (y,s) R x [-Mi, Mi] and 
< h < 1 we have 

f(Y h , S h ) - /(y, s) 



h 



0. 



We may, therefore, assume that (y,s) G R x [—Mi, Mi]. We will also assume that 
< h < 1. 

As above we see that for any starting point (Yq, So) = (y, s) G R x [—Mi, Mi] and all 
< t < 1 we have \S t \ < M x + WW^. Set M 2 := M 1 + ||W^"||oo- We assume without 
loss of generality that X = 0. Then 



Y t = y + X t + [ W'(Y r )S r dr, 
Jo 

S t = s+ [ W"(Y r )dr. 
Jo 



It follows that 

T h f(y, s) - f(y, s) E^) f(Y h , S h ) - f(y, s) 



h 



h 



l - E^[f(Y h , S h ) - f(Y h , s)} + l - E^[f(Y h , s) - f(y, s)} 



= 1 + 11. 

Using Taylor's theorem we find 



I = E (j/,s) 



19/ 

h ds 



h x 2- 

W"(y r ) dr 



l% {Yh > S) Jo W " {V) dT + vfs {Yh ' S) fo {W " {Yr) ~ W " {V)) dT 
(Y h ,0 (^£w"(Y r )dr 



2hds 2 
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where £ is a point between s and Sh- Note that 

-h 



l^ {Yh,s) l ( w "( Y r)- w "(y)) dr 



df 



df 



< 



ds 
*0, 



ds 



X r + / W'(Y t )S t dt 
Jo 



A2\\W"\\n > dr 



sup \X r \ + /i||W||ooM 2 ) } A 2||W"| 

0<r</i 



uniformly for all (y,s) G R x [—Mi, Mi]. The convergence follows from the right 
continuity of X t and our assumption that Xo = 0. We also have 

2 



1 d 2 / 



2/i 9s 2 



(r fc ,0 (£w"(Y r )dr 



< 



d 2 f 



ds 2 



-\\W"\\ 2 

2 11 "°° 



h->0+ 



+ 0, 



uniformly for all (y,s) G R x [—Mi, Mi]. Because Y h is right-continuous it is easy to 
see that 



if (W W» (v)i r 



^(y,s)W"(y), 



h->0+ ds 

uniformly for all (y, s) G R x [-Mi, Mi]. It follows that 

uniformly for all (y, s) G R x [—Mi, Mi]. 
Now let us consider II. We have 

II = I E^[f(y + X h , s) - f(y, s)] + i E<**>[/(r fc , s) - f(y + X,, s)] 

= Hi + H 2 . 
It is well known that 

ih-—>-(-\T /2 f(y,s), 



uniformly for all (y, s). We also have 

'ldf rh 
h dy 



II, = E (y ' s) 



(y + X h ,s)J W'(y r )S r dr + ^0(£,s)Qf W"(y)S r dr 

(y + x,, s ) ( jT ( y ) s rfr + jT" ^'(y)(^ - s) dr 

+ jf (W'(y r ) - W\y))sdr 



13/ 



1 «9 2 f 



(£, s ) / vy'(y)5 r dr 



where £ is a point between y + X h and Y^. Using similar arguments as above we obtain 

H 2 



>^(y,a)W / (y) a , 
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uniformly for all (y, s) G R x [—Mi, Mi]. 
It follows that 

T h f(y,s)-f(y,s) -(-A^fiv, s) + W\y)s 9 My, s) + W"(y) & My, s), 

h h->o+ ay as 

uniformly for all (y,s) G R x [-Mi, Mi]. This means that / G r D(9 (Y ' s) ) and (j277jl 
holds. □ 



Remark 2.6. A weaker version of Lemma 12.51 can be proved as follows. If we rewrite 
the SDE (J2ZD in the form 

t 



<g)-(SW)'(?) -•«■*>' 



and notice that (X 4 ,t) T is a two-dimensional Levy process with characteristic expo- 
nent ip(£,, T ) — \£\ a + ? " T > we can use [H3 Theorem 3.5, Remark 3.6] to deduce that 
C^°(R 2 ) C D(S (T ' S ' ) )- This argument uses the fact that the SDE has only jumps in 
the direction of the a-stable process, while it is local in the other direction. The- 
orem 3.1 of [15] now applies and shows that S^ ¥ ' s ^ is a pseudo-differential operator 
2 Y,s u(x, s) = (27r) -2 J r2 p(x, s; £, r) ?«(£, r) e lx ^ +tST dr, where 5" denotes the Fourier 
transform, with symbol 

p(x, s; £, r) = <K<%, s) T (£, r) T ) = + i£W'(x)s. 

A Fourier inversion argument now shows that (12. 7p holds for / G C^°(R 2 ) and by a 
standard closure argument we deduce from this that (12. 7p also holds for / G Cq(R 2 ). 

We say that / G C (S x R) if and only if for every e > there exists a compact set 
K C S x R such that \f(u)\ < e for u G K c . Let us define the semigroup {Tf} t >o of 
the process (Z t , S t ) by 

(2.8) T?f(z,s) = E^f(Z t ,S t ), zeS, s G R, 

for / belonging to Cq(S x R). Let z = e iy , y G R. For future reference, we note the 
following consequences of Corollary 12. 3\ 

(2.9) 7f/(z, s) = E^ f(Z t , S t ) = E^ f(j Y *,S t ) = E^ f(Y t , S t ) = T t f(y, s), 
and 

(2.10) Tff(y,s) = Tj(y,s). 

By Arg(z) we denote the argument of z G C contained in (— n, 7r]. For g G C 2 (S) let 
us put 

Lgfz) = A a lim / f ^-j — dio 

^ 0+ isn{|Arg(W^)l> e } |Arg(w/^)| 1 + Q 



(2.11) 

d \ / - 

Arg(w/z) + 2mr| 1+Q 



neZ\{0} 1 



where .Aq, is the constant appearing in ( 12. 5ft and c?w denotes the arc length measure on 
S; note that J s (iu; = 2n. 
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Let S be the generator of the semigroup {Tf }i> and let 2)(S) be its domain. 
Lemma 2.7. We have C 2 ($ x E) C D(9) and for f G C 2 ($ x R), 

Sf(z, s) = L z f(z, s) + V'(z)sf z (z, s) + V"(z)f s (z, s), z G S, sGi 

Proo/. Let / G C 2 (S x E). Note that / G C 2 (E 2 ). We obtain from (EZJ), for 2 = e iy , 
y, s g R, 

nm 1*f(z,a)-f(z,s) = Um Jj /(y, g) - f(y, a) 
t^o+ t t->o+ t 

(2-12) = -(-A)«/ 2 /(y, a) + s) + ^"(y)/ s (y, a). 

By Lemma [2.51 this limit exists uniformly in z and s, i.e. / G D(S). 
We get from (12.51) 

(2.13) - {-\T /2 f(y,8) = L,f(z,8). 

Recall that we have W(y) = V(e iy ), y G E. Using our definitions we get V'(z) = W'(y), 
V"(z) = W"{y) for z = e iy , y G E. Hence fl2TT2l) equals 

L z f(z, s) + V'(z)sf z (z, s) + V"(z)f s (z, s), 

which gives the assertion of the lemma. □ 

We will need the following auxiliary lemma. 

Lemma 2.8. For any f G C 2 (S) we have 

Lf(z) alz = 0. 



J$ 

Proof. Recall that Arg(z) denotes the argument of z G C belonging to (—71,71]. First 
we will show that 

(2.14) J 1 {W : \Avg(w/z)\>e}{w) ^^J^ dw dz = 0. 

We interchange z and w, use Fubini's theorem and observe that | Aig(z/w) \ = \ Arg(w / z) 
LJ {W : ' A ^ w/Z)l>e}{w) \Arg(w/z)\^ dW dZ 

[[ t m-f(w) 




SxS 

f 

SxS 



which proves (12.140 . 




I { 2 : | Arg(z/«,) | >e} [?) J^JJ^j 1 1 +Q dZ dW 

f(z)-f(w) 

I{ Z : | Arg(^)|> £ }U) | Arg(2/w)| 1+Q 

, x f{w)~f(z) , , 
1{ W : | ArgH/ 2 )|> £} H |Arg(w/z)|1+Q ^ <fc, 
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By interchanging z and w we also get that 

(2 ' 15) = y II m^m dzdw 



Note that for Arg(w/2;) ^ n we have | Arg(z/w) + 2nir\ = \ Arg(w/z) — 2mr\. Hence 
the expression in (12.151) equals 0. 
Set 

Lsf(z):=l dw. 

Jsn{\Arg(w/z)\>e} \ Alg{w / z)\ 1+a 

What is left is to show that 



(2.16) / lim L £ f(z)dz = lim LJ(z)dz. 

JS e-»-0 + Js 

By the Taylor expansion we have for / G C 2 (S) 

/(«;) - /(*) = Axg(w/z)f'(z) + Arg 2 (w/2;)r(w, z), w.zGS, 
where |r(iu, ,z)| < c(f). Hence, 

\L £ f(z)\= I r{w,z)kvg 1 - a {w/z)dw 

JsC\{\Arg{w/z)\>e} 

<c(f) I \Arg 1 - a (w/z)\dw = c(f,a). 
Js 

Therefore, we get H 2 . 1 6 1) by the bounded convergence theorem. □ 

We will identify the stationary measure for (Z t , S t ). 
Proposition 2.9. For z G S and s G R let 

1 1 2/ 

p 1 (z) = — , p 2 (s) = — p= e" s /2 , ir(dz,ds) = p 1 (z)p 2 (s) dzds. 

27T V27T 



Tnen /or any / 6 C C 2 (S x R) we have 

/ / 9f(z,s)Tc(dz 1 ds) = 0. 

JS JR 

Proof. We have 

/ / 9f(z, s)ir(dz,ds) 
Js Jr 

= -!-/ I (L z f(z,s) + V'{z)sf z {z,s) + V"(z)f s (z,s))p 2 (s)dsdz. 
Z7r Js Jr. 
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Integrating by parts, we see that this is equal to 

— [ L z f(z,s)p 2 (s)dsdz - — / / V"(z)sf(z, s)p 2 {s)dsdz 
27T Js Jr 271 Js Jr 

~ [ f V"{z)f{z,s)p 2 '{s)dsdz = 1 + 11 + III. 
27r Js Jr 

Since p 2 '(s) = —sp2{s) we find that II + III = 0, while I = by Lemma [2.81 The claim 
follows. □ 



Proposition 2.10. For any t > we have 

T? : C 2 (S x R) ->■ C7 C 2 (S x R). 

The proof of this proposition is quite difficult. It is deferred to the next section in 
which we prove this result in much greater generality for solutions of SDEs driven by 
Levy processes. 

Theorem 2.11. Let 

(2.17) ir(dz,ds) = . * e~ s2/2 dzds, z G $, s G R. 

(2n) i ' 2 

Then ir is a stationary distribution of the process (Z t ,S t ). 

Proof. Let (Y t , S t ) be a Markov process satisfying the SDE ( 12. ip and (Z t , S t ) = (e lYt , St). 
Recall that {T®} t > is the semigroup on C (S x R) defined by (I2.8P and S is its generator. 
Let CP(R x R) and CP(S x R) denote the sets of all probability measures on R x R and 
S x R respectively. In this proof, for any p G ^(S x R) we define p G T(R x R) by 
p{[0, 2tt) x R) = 1 and p{A x B) = p{e iA x B) for Borel sets A C [0, 2vr), B C R. 

Consider any p G CP(S x R) and the corresponding p G CP(R x R). 

For this p there exists a Markov process (Y t , S t ) given by (12. ip such that (Y , So) has 
the distribution p. It follows that for any p G CP(S x R) there exists a Markov process 
(Z t , S t ) given by (12. ip and Z t = e tYt such that (Z Q , S ) has the distribution p. By 
Proposition 4.1.7 [TU], (Z t , S t ) is a solution of the martingale problem for (S,pO- The 
Hille-Yosida theorem shows that the assumptions of Theorem 4.4.1 [10J are satisfied if 
we take A = A' = 9- Thus Theorem 4.4.1 [TU] implies that for any p G CP(S x R), 
uniqueness holds for the martingale problem for (S,p). Hence the martingale problem 
for S is well posed. 

Note that C 2 (S x R) is dense in Co(S x R), that is, in the set on which the semigroup 
{Tf}i> is defined. It follows from Proposition 12. 101 and Proposition 1.3.3 from [10] that 
C 2 (S x R) is a core for 9- Now using Proposition 12.91 and Proposition 4.9.2 from [10] 
we get that ir is a stationary measure for 9- This means that (Z t , S t ) has a stationary 
distribution ir. □ 



Theorem 2.12. The measure ir defined in (I2.17P is the unique stationary distribution 
of the process (Z t , S t ). 
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Proof. Suppose that for some cadlag processes X 1 and X 2 , processes (Y^Sj) and 
(Y 2 , S 2 ) satisfy 

(2.18) Y t 1 =y + X] + f W'(Y r l )Sl dr, 

Jo 

(2.19) Si = 8+ I W{Y})dr, 

Jo 

(2.20) Y t 2 = y + X 2 + f W'(Y 2 )S 2 dr, 

Jo 

(2.21) S 2 = s+ [ W"{Y 2 )dr. 

Jo 

Then 

(2.22) \S} - S 2 \ < f \W"{Y?) - W"{Y 2 )\ dr < \\W^\U f \Y} - Y 2 \ dr, 

Jo Jo 

and, therefore, for t < 1, 

\Y t ' - Y 2 \ < \Xl - X 2 \ + f \W'(Y})Sl - W\Y 2 )S 2 ] dr 

Jo 

< \Xl - X 2 \ + f \W'{Y^){Sl - S 2 )\ dr + f \(W'tf) - W'(Y 2 ))S 2 \ dr 

Jo Jo 

< \Xl-X 2 \ + \\W'\\ 00 snp ^l-S^t+WW'W^ sup \S 2 r \ C \Yl-Y 2 \dr 

0<r<t 0<r<t Jo 

< \Xl - X 2 \ + UWOof Halloo f \Y r l - Y 2 \ dr 

Jo 

+ \\W"\\ 00 (\s\ + \\W"\\ 00 t) f\Y^-Y 2 \dr 

Jo 

< \xl - x 2 \ + ( Cl t + c 2 \s\) f\Y}-r; 

Jo 



| dr 
dr. 



<\Xl-X 2 \ + { Cl + c 2 \s\) [ \Y r l -Y 2 

Jo 

By Gronwall's inequality, 

sup \Y^-Y 2 \< sup \X^-X 2 \ + / \X l r -X 2 \{ Cl + c 2 |s|)exp{(ci + c 2 \s\)t}dr 

0<r<t 0<r<t Jo 



< sup iX^-X^Kl + ^ci + CalsDexpl^i + Calsl)*}). 

0<r<t 

For t — 1, the inequality becomes 

(2.23) sup \Y* - Y 2 \ < sup \X l r - X 2 r \(l + (ci + c 2 |s|) exp {(a + c 2 |s|)}). 

0<r<l 0<r<l 
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We substitute (12.191) into (12.181) and rearrange terms to obtain, 

X] = -y + Yl - jf* W'O?) (a + jf W"(Y^) du^j dr. 
We substitute the (non-random) number y for in the above formula to obtain 
(2.24) xl = -y + y- jf W\y) + jf W"^) du) dr 

= -^'(y)(t S + t 2 iy"(y)/2). 

From now on, X 1 will denote the process defined in (12.241) . It is easy to see that X\ is 
well defined for all t > 0. If we substitute this X 1 into (l2TT8l) - (T2~T9l ) then K t = y. 

It follows from [T8J Theorem II, p. 9], that every continuous function is in the support 
of the distribution of the symmetric a-stable Levy process on R. We will briefly outline 
how to derive the last claim from the much more general result in [LSI Theorem II, p. 9]. 
One should take a( ■ ) = and b( ■ , z) = z. Note that the "skeleton" functions in [T8J (5), 
p. 9] can have jumps at any times and of any sizes so the closure of the collection of all 
such functions in the Skorokhod topology contains the set of all continuous functions. 
Standard arguments then show that every continuous function is in the support of 
the distribution of the stable process also in the topology of uniform convergence on 
compact time intervals. We see that if X 1 is the continuous function defined in ( I2.24p 
and X 2 is a stable process as in ( 12. ip then for every e > there exists 5 > such that, 



P ( sup \X l r 

\0<r<l 



X 2 r \ < ej > 5. 

This and (I2.23P show that for any y, s € R and s > there exists 5 > such that, 



( sup \X] 


— x 2 | 


< e, sup F r 2 — y\ 


\0<r<l 




0<r<l 



Note that S can change by at most HW'Hoo on any interval of length 1. This, the 
Markov property and induction show that for any e > there exist 6% > 0, k > 1, such 
that, 



( sup \X] - X*\ < 2- k e, sup \Y r 2 - Y*\ < 2- k e] > 5 k . 

\fc<r<fc+l fc<r<fc+l / 



where X 1 is defined in ( I2.24p . This implies that for any r < oo, y, s £ R and e > 
there exists 5' > such that, 

(2.25) F y ' s ( sup \Xl - X 2 | < 2e, sup \Y r 2 - y\ < 2e j > 5'. 

\0<r<T 0<r<T J 

Step 2. Recall that V is not identically constant. This and the fact that V £ C 5 easily 
imply that W" is strictly positive on some interval and it is strictly negative on some 
other interval. We fix some ai,a,2 £ (—n,n), b\ > 0, b 2 < and Eq £ (0,7r/100), such 
that V"(z) > bi for z £ S, Arg(z) £ [a x - 4e ,ai + 4e ], and V"(z) < b 2 for z £ S, 
Arg(^) £ [a 2 - 4e , a 2 + 4e ]. 
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Suppose that there exist two stationary probability distributions tt and 7? for (Z, S). 
Let ((Z t , S t ))t>o and ((Z t , S t ))t>o be processes with (Z , So) and (Z , So) distributed 
according to ir and 7?, respectively. The transition probabilities for these processes are 
the same as for the processes defined by ( 12. ip and ( 12. 4ft . Let X denote the driving 
stable Levy process for Z. 

Let A be an open set such that W"(y) > c > for all y G A. In view of the 
relationship between V and W, we can assume that A is periodic, that is, y E A if and 
only if y + 27r G A. It follows easily from (12. ip that there exist q% > and si < 00 
such that for any (Y , So), the process y enters A at some random time T\ < s\ with 
probability greater than q 1 . Since Y is right continuous, if Y Tl G A then y t stays in 
A for all £ in some interval (Ti,T 2 ), with T 2 < 2s x . Then (12. ip implies that S t ^ 
for some £ G (Ti,T2). A repeated application of the Markov property at the times 
2si,4si,6si, . . . shows that the probability that St — for all £ < 2ks\ is less than 
(1 — qi) h . Letting k — » 00, we see that St ^ for some £ > 0, a.s. 

Suppose without loss of generality that there exist e\ > 0, £ 2 > and pi > such 
that P 7r (S' t2 > £1) > pi. Let F 1 = {S t2 > £1} and t 3 = £i/(2||iy"|| 00 ). It is easy to see 
that for some p 2 > 0, 

(3 1 G [t a , t 2 + £3] : Arg(Z t ) G [a 2 - £ 0) a 2 + e ] | Fj > p 2 . 

This implies that there exist £1 > 0, t 2 > 0, t 4 G [t 2 , t 2 + £3] and p 3 > such that, 

F n (S t2 > £1, Arg(Z t4 ) G [a 2 - 2£ , a 2 + 2e )) > p 3 . 

Note that \S U - S t2 \ < ||W"||oo*3 < £i/2. Hence, 

P^ 4 > £l /2, Arg(Z t4 ) G [a 2 - 2£ , a 2 + 2e ]) > p 3 - 

Let £ 2 G (£i/2, 00) be such that 

F*{S U G [£i/2,£ 2 ], Arg(Z t4 ) G [a 2 - 2£ ,a 2 + 2e }) > p 3 /2. 

Let £5 = 2£ 2 /|6 2 | and ts = £4 + £5. By (I2.25p . for any £3 > and some p± > 0, 

( sup \X} - A r | < £ 3 , S u G [£i/2, £ 2 ], 

Aig(Z t ) G [a 2 - 3£ , a 2 + 3£ ] for all £ G [£ 4 , U] \ > p 4 , 

where X 1 is the function defined in ( I2.24p . Since V"(z) < b 2 < for Arg2 G [a 2 — 
3eq, a 2 + 3£o], if the event in the last formula holds then 

S t6 =S t4 + ^ V"(Z S ) ds<e 2 + b 2 t 5 < -£ 2 . 

This implies that, 

(2.26) P*( sup \X] - X r \< £ 3 , S t4 > £i/2, Ste < -£2) > Pa- 

^t 4 <r<t 6 ' 

Step 3. By the Levy-Ito representation we can write the stable Levy process X in the 
form Xt = Jt + X t , where J is a compound Poisson process comprising all jumps of X 
which are greater than £ and X = X — J is an independent Levy process (accounting 
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for all small jumps of X). Let us denote by A = A(cu, £o) the rate of the compound 
Poisson process J. 

Let (Y, S) be the solution to ( 12. ip . with X t replaced by X t for t > t 4 . Take £ 3 < e /2. 
Then sup i4<r<if . \X} — X r \ < £ 3 entails that sup t4<r<t(3 \J ti — J r \ = 0. Thus, ( I2.26P 
becomes 

P w ( sup \X l r - X r \ < £ 3 , S u > f, S t6 < -e 2 ) 

>P"( sup \X}-X r \<e 3 , sup \J u -J r \ = 0,S t4 >f,S t6 <-e 2 ) 

M 4 <r<i 6 U<r<t 6 ' 

> pi > 0. 

Let t be the time of the first jump of J in the interval [t4, tg] ; we set r = t§ if 
there is no such jump. We can represent {(Y t ,S t ),0 < t < r} in the following way, 
{Y t , S t ) = {% S t ) for < t < r, S T = S T , and Y T = Y T + J T - J T _. 

We say that a non-negative measure /ii is a component of a non-negative measure \x 2 
if \i 2 = for some non-negative measure ^3. Let fi(dz, ds) = F^^Zr G dz, S T G ds). 

We will argue that fi(dz, ds) has a component with a density bounded below by c 2 > 
on S x (—£2,61/2). We find for every Borel set A C S of arc length \A\ and every 
interval (si,s 2 ) C (— £ 2 ,£i/2) 

fi(A x (si,s 2 )) 
= P^(Z r GA S T e(si,s 2 )) 

> (z T eA,S T e (s u 82 ), sup \X l r -X r \< £ 3 , 5 i4 > f , S te < -£ 2 

^ ti<r<te 



> F" (V (JW ^ } G e'^'A, S T G (si,s 2 ), 



sup \Xl -X r \< £ 3 , 5 i4 > £x/2, S t6 < -£ 2 , N J = 1 

i4<r<t6 

Here iV J counts the number of jumps of the process J occurring during the interval 
[^4,^6]- Without loss of generality we can assume that £0 < 2n. In this case the density 
of the jump measure of J is bounded below by c 3 > on (2n , 4n) . Observe that the 
processes (X,S) and J are independent. Conditional on {iV J = 1}, r is uniformly 
distributed on [t 4 , t 6 ], and the probability of the event {N J = 1} is A(t 6 — t 4 )e" A( * 6 "* 4) . 
Thus, 

fi{A x (si,s 2 )) 

>c 3 |A| (S T G ( Sl , s 2 ) sup \X l r -X r \< £3, S u > £l/2, S t6 < -£ 2 , N J = l) 
x p 4 ■ A(t 6 -t 4 )e" A( * 6 "' 4) . 

Since the process S spends at least (s 2 — si)/||W"||oo units of time in (s 1; s 2 ) we finally 
arrive at 

(si, s 2 )) > p,Xe- x{te - u) c 3 \A\(s 2 - sO/IIW'lU. 
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This proves that fi(dz, ds) has a component with a density bounded below by C2 = 
p 4 A e -^-*4) C3 /||W|| 00 onSx(-e 2 ,e 1 /2). 

Step 4- Let £4 = £i/2 A £2 > 0. We have shown that for some stopping time r, 
P 7r (Z T G dz,S T G <is) has a component with a density bounded below by C2 > on 
S x (—£4, £4). We can prove in an analogous way that for some stopping time r and 
£4 > 0, F n (Zr G dz, St G ds) has a component with a density bounded below by 02 > 
on S x (—£4, £4). 

Since tt 7^ 7?, there exists a Borel set A C S x R such that vr(y4) 7^ tt(^4)- Moreover, 
since any two stationary probability measures are either mutually singular or identical, 
cf. [HI Chapter 2, Theorem 4], we have tt(A) > and tt(A) = for some A. By the 
strong Markov property applied at r and the ergodic theorem, see [191 Chapter 1, page 
12], we have P^-a.s. 

lim (1/t) / l {( z s ,5 s)e A } ds = ir{A) > 0. 
Similarly, we see that P^-a.s. 

lun(l/t)^l {( ^ §3)eA} rf S = vr(A) = 0. 

Since the distributions of (Z T , S T ) and (Zf, SV) have mutually absolutely continuous 
components, the last two statements contradict each other. This shows that we must 
have tt = tt. □ 

Remark 2.13. It is not hard to show that Theorem 12.111 holds even if we take a = 2 
in ( 12. ip . that is, if X t is Brownian motion. It seems that for a = 2 uniqueness of the 
stationary distribution can be proved using techniques employed in Proposition 4.8 in 
[TJ. A close inspection of the proofs in this section reveals that our results remain also 
valid if X t is a symmetric Levy process with jump measure having full support. 



3. Smoothness of T t f 

In this section, we will show that if / G Cf then T t f G where {T t }t>o is the semi- 
group of a process defined by a stochastic differential equation driven by a Levy process. 
We use this result to show Proposition 12. 101 but it may well be of independent interest. 
We found some related results in the literature but none of them was sufficiently strong 
for our purposes. The key element of the proof are explicit bounds for derivatives of 
the flow of solutions to the SDE. This is done in Proposition 13.31 We provide a direct 
and elementary proof of this proposition. Note that our bounds are non-random and 
do not depend on the sample path. This is a new feature in this type of analysis since 
usually, see e.g. Kunita [T2], the constants are random since they are derived with the 
Kolmogorov-Chentsov-Totoki lemma or a Borel-Cantelli argument. Let us, however, 
point out that there is an alternative way of proving Proposition 13.31 It is possible 
to use [IU Theorems V.39, V.40] and [HJ formula (D), p. 305] to obtain bounds for 
derivatives of the flow. Since this alternative approach demands similar arguments and 



18 K. BURDZY T. KULCZYCKI AND R.L. SCHILLING 

is not shorter than our proof of Proposition 13.31 we decided to prove Proposition 13.31 
directly. 

Consider the following system of stochastic differential equations in R n , 

' dY x {t) = dX x (t) + VxiYitydt, 

(3.i) | ; 

dY n (t) = dX n (t) + V n (Y(t))dt, 

where Y(t) = (Y^t), . . . ,Y n (t)) E R n , X(t) = (X^t) , . . . , X n (t)) E R n . We assume 
that X(0) = 0, Xi, . . . ,X n are Levy processes on R and Vi : R n — > R are locally 
Lipschitz. We allow Xi, . . . , X n to be degenerate, i.e. some or all Xi may be identically 
equal to 0. 

By [HI Theorem V.38] it follows that if Y(0) = x then there exists a stopping time 
((x, uj) : R n x Q — )■ [0, oo] and there exists a unique solution of (13.1 ft with Y(0) = x with 
lim sup^^.) \ Y(t)\ = oo a.s. on ( < oo; ( is called the explosion time. In order to apply 
[HJ Theorem V.38] we take in the equations marked (®) in [H,, p. 302] m — n + 1, 
Xi = Yi(t), x* = Yi(0), Zf = X a (t) for a E {1, . . . , n}, Z? +1 = t and /» = 5 ai for 
a, i E {1, . . . ,n} and f % n+l (x) = Vi(x) for i E {1, . . . , n}. 

By Y x (t) we denote the process with starting point Y x (0) = x. In the rest of this 
section, we will assume that ( 13. II) holds not only a.s. but for all u E Q. More precisely, 
we can and will assume that the solution to ( 13. ip is constructed on a probability space 
Q such that X(0) = and 

Y x (t) = x + X(t)+ [ V{Y{s))ds, 



for all £ > and all u E Q. 

Set 

|| x || — max{ \ x± | , . . . , \x n | } , x — {x\ , . . . , x n ) , 

and 

B*(x,r) = {y E R" : \\y - x\\ < r}, x E R n , r > 0. 
For / : R n -> R and A C R n we write L>W/ = V/, 

H/IUa = sup |/(a;)|, pW/lloM = E SU P l^/COl, 

||/|| -),A = ||/||oo,A + P (1) /l|oo,A + • • • + P (i) /l|oc,A. 

When A = R n we drop A from this notation. For V = (Vi, . . . , V n ) from (13. ip and 
A C R n we put 

n n 
i=l i=l 

IMIo-),a = ll^iu^ + IP^Vlu,^ + . . . + \\d^v\u a . 

For / : R™ — > R, x E R n and < t < oo we define the operator T t by 
(3.2) T t f(x)=E[f(Y*(t));t<((x) 
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Before formulating the results for the process Y(t) let us go back for a moment to 
the original problem (12. ip . that is, 

'dY t = dx t + w'(y t )s t dt, 

dS t = W"(Y t )dt. 
This SDE is of type (13. ip because we can rewrite it as 

'dYi(t) = dXi(t) + V[(y(t)) dt, 
dY 2 {t) = dX 2 {t) + V 2 (Y(t)) dt, 



(3.3) 



where Xi(t) = X t is a symmetric a-stable Levy process on R, a G (0,2), X 2 {t) = 0, 
V\ (2/1,2/2) = W'(yi)y 2 , V 2 (yi,y 2 ) = W"(yx). By Lemma IO there exists a unique 
solution to this SDE and the explosion time for this process is infinite a.s. We want 
to show that T t f G C 2 whenever / G C 2 . Our proof of Theorem 13.11 requires that V{ 
and its derivatives up to order 3 are bounded. However, V\ (2/1,2/2) = W(yi)y 2 is not 
bounded on R 2 . We will circumvent this difficulty by proving in Proposition 13.61 that 
T t f G C 2 (R 2 ) whenever / G C 2 (R 2 ), where C*(R 2 ) is given by Definition E2J 

Let us briefly discuss the reasons that made us choose this particular set of functions, 
C 2 (R 2 ). This discussion gives also an explanation for the specific assumptions in the 
main result of this section, Theorem 13.11 

Assume that / G C 2 (R 2 ) and supp / C K = R x [—r,r], r > 0. Fix t < 00. If 
\s\ = \S \ > r + toll W'lloo then for t < t , 



s + 



W"{Y^)du 
Ef(Y t ^,S^) = 



> r 



0. 



3 , r + to||W"||oo]- 



and, therefore, 

T t f(y,s) 

It follows that if t < to then 

(3.4) supp(T t /) C K = R x [-r-t \\W"\^. 

For technical reasons, we enlarge K as follows, 

K 3 — Rx ( — r — to||W"||oo - 3, r + t^W"^ + 3). 

In view of (13. 4p . we have to consider only starting points (y, s) G K in order to prove 
that T t f G C 2 (R 2 ). Note that for the starting point (y, s) G K 3 and t <t Q we have 



SI 



W"{Y^ s) )du 



<r + 2t ||^"| 



3. 



Thus for all starting points (y, s) G K 3 and t < to, 

(3.5) (Y t iv ' s) , Sl y ' s) ) GM:=Rx [ - r - 2* ||W /, || O a - 3, r + 2t \\W"\\ oa + 3] . 

But the function Vx(yi,y 2 ) = W'(yi)y 2 is bounded on M. Using our assumptions on 
W, namely, periodicity of W and 1^ 6 C 5 , we obtain also that the derivatives of 
V\ (2/1,2/2) = W{i/i)y 2 up to order 3 are bounded on M. 
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Now we return to the general process Y(t). Let us formulate the main result for this 
process. 

Theorem 3.1. Let f : R n ->■ R be a function in C$. Fix < t < oo. Let Y x (t) 
be a solution of (13. ip . Assume that the explosion time ((x,u) = oo for all x G R n 
and all u G Q. Let T t f be defined by (13.21) . Assume that K C R n , for every i < in 
supp(T 4 /) C K and that there exists a convex set M C R n such that Y x (t,u>) G M 
for all x G K% := [j x&K B*(x,3) , t < t , and u G Q. Assume that HV^I^m < oo and 
\\D^V\\ao,M < 00 f or 3 — 1,2,3. Then we have 

T t f G C 2 h for all t < t . 

Remark 3.2. When ||V||(3) < oo (i.e. when the assumptions of Theorem 13.11 hold with 
K = M = R") then the above theorem implies that we have for any / G 

T t feCl for all t > 0. 

The first step in proving Theorem 13.11 will be the following proposition. 

Proposition 3.3. Fix < to < oo. Let Y x (t) be a solution of (13.11) . Assume that 
the explosion time ((x,u) = oo for all x G R n and all u G Q. Let K C R n . Assume 
that there exists a convex set M C R n such that Y x (t,u) G M for all x G K 3 := 
[J K B*(x,3), t < t , and u G Q. Assume that WVW^^m < oo. Put 

/ \ 1 /l 

(3.6) r := — - — 7- — A in, - := oo 

For every u G £1 we have the following. 

(i) For allO <t <r, x G K 2 = \J x€K B*(x, 2), h G R n , \\h\\ < 1, 

(3.7) \\Y* +h (t,u)-Y*(t,u)\\ <2\\h\\. 

(ii) Recall that ej is the i-th unit vector in the usual orthonormal basis for R n . For 
all < t < r , x G K2, i G {1, . . . , n}, 

exists, and 

(3.8) \\D l Y x (t,uj)\\<2. 

We will write DiY x {t, u) = (A^(i, w), . . . , DiY x (t, u)). 
(hi) For allO <t <t, x G K x = \J xeK B*(x, I), h G R n ; \\h\\ < 1, % G {1, . . . , n}, 

(3.9) \\DiY x+h {t,u) - DiY*(t,u)\\ < 8 WD^VW^m r \\h\\. 
(iv) For all < i < r, x G K x , i, k G {1, . . . , n}, 

D ik Y (t,u) := lrm 

u^O U 

exists and 

(3.10) \\D ik Y x (t,u)\\ ^SWD^VW^mt. 

We will write D ik Y x (t,u) = (D lk Y x (t,u), . . . , D ik Y*{t,uj)). 
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(v) For all < t < r, x G K , h G R n , \\h\\ < 1, i, k G {1, . . . , n}, 
\\D lk Y x+h (t,co)-D lk Y x (t,u)\\ 

< 96 WD^VWl^r 2 \\h\\ + 16 {{D^V^m t \\h\\. 

Remark 3.4. The existence of DjY x (t) and D ik Y x (t) follows from [TJl Theorem V.40]. 
What is new here are the explicit bounds for DiY x (t) and D ik Y x (t) which are needed in 
the proof of Theorem 13.11 see Lemma 1331 The proof of Proposition 13 .31 is self-contained. 
We do not use [HI Theorem V.40]. 

Proof of Proposition Iff. 31 The proof has a structure that might be amenable to presen- 
tation clS 8b CclSG of mathematical induction. After careful consideration we came to the 
conclusion that setting up an inductive argument would not shorten the proof. 

Recall that we assume that ( 13.1 j) holds for all u G ft, not only a.s. Throughout this 
proof we fix one path u G £1. 

(i) Let x G K 2 , h G R n , \\h\\ < 1 and < t < r. Recall that X(0) = 0. For any 
1 < j <nwe have 

(3.11) Y x+h (t) - Y x (t) = hj+ [ [V 3 {Y x+h {s)) - V 3 {Y x {s))] ds. 

Jo 

Let 

Cl := Cl (x,h) := sup \\Y x+h (t)-Y x (t)\\. 

0<t<T 

Note that for < t < r we have Y x (t) G M and Y x+h (t) G M. By (13TTT]) and 
||V||oo,m < oo we get that C\ is finite. Moreover, 

\\Y» +h (t)-Y x (t)\\ < \\h\\+ [ t \\D^V 3 \\ ooM \\Y x+h (s)-Y x (s)\\ds 

Jo 

< \\h\\ + tWD^VjW^mCl 

Hence, 

ci < \\h\\ + t||L> (1) V|| 00 , m c 1 , 
which, when combined with ( I3.6p . gives 

sup \\Y x+h (t) - Y x (t)\\ = Cl < - J^L < 2\\h\\. 



(ii) Denote 

"3 \"J 3 \ ) j 

-,x,h 



R x > h (t) = Y x+h (t)-Y x (t) 



and R x ' h {t) = (R x ' (£), . . . , R x ' h (t)). Using the Taylor expansion we get from (13. lip . 

(3.12) Rf h (t) = h j + [ D^V j (Y x (s))-R x ' h (s)ds + 0(\\h\\ 2 ). 

Jo 

For i G {1, . . . , n} and h = ue^ let 

/ s /V *?*(*) ,. • ,. «?*(*)' 

c 2 = c 2 (z,z) = max sup hmsup — hmmf — 

!<?'<" 0<t<r \ u^O U u->0 u 
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Note that c 2 is finite because for u G (—1, 1) we have \R x,h (t)\ < 2u, by (13. 7L Consider 
< t < t, x G K 2 , i,j G {1, . . • ,n}. From (I3.12p we obtain for u,u' G (-1, 1) \ {0}, 
h = uti and h' = u'ei, 

Rf(t) Rf(t) _ f » ( s ) tff ( s ) 



K It' 

Letting it, w' — >• leads to 

limsup -^-ll - liminf J , < r L> (1 V U,M • c 2 , 

and since < t < t and j G {1, . . . , n} are arbitrary, we get 

c 2 <t\\D^V\\ OC}M -c 2 . 

So C2 = which means that DiY x (t) exists. Estimate (I3.8P is now an easy consequence 
of CT- 

(iii) From f 1 3 . X 2 j) and the bounded convergence theorem, we obtain 

(3.13) DiY?(t) = 5^ + f D^V j (Y x (s)) ■ D t Y x (s)ds. 

Jo 

Let x G K u he R n , \\h\\ < 1 and i G {l,...,n}. Set 

c 3 :=C3(z,M):= sup || DiY x+ \t) - A^(t)||- 

0<t<T 

Because of (13.81) . C3 is finite. For any < t < r we have 
DiYf+^t) - DiY?(t) 

[ [D^VjiY^is)) ■ DiY x+h (s) - D^Vi(Y x (s)) ■ DiY x (s)] ds 
Jo 

[D^VjiY^is)) - D^VjiY^s))] ■ DiY^s) 
+ D^V 3 (Y x (s)) • [Af +ft (s) - DX x (s)] )ds, 



(3.14) 



so 



|A>7 +A (t)- A*7(*)| < / 

Jo 



E 1 A^(r^(,)) - ^r( S ))i lA^r^^) 



fc=i 

n 



+ y £\D k V j (Y*(s))\ \DiY x+h (s) - DiY£(s) 

k=l 

In view of (13.71) and (I3.8p . we have for < s < r, 

n 

^2\D k V j (X^ h (8))-D k Y l (Y'(8))\ < \\D^V\\ ooM \\y X+h is)-Y x (s)\\ 

k=l 

<2|| J D( 2 )V|| 0O , M | 



ds. 
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n 

IIAl^OOII < 2, \DkVj(Y*(s))\ < WD^VUm. 

fc=i 

It follows that 

\DiYj? +h (t) - DiYj°(t)\ < 4 WD^VW^m r \\h\\ + r H^VIUm • c 3 , 

so, 

oo,M T \\h\\ + T ||-D^y||oo,M " C3. 

By definition, r < l/(2p( 1 )V r || 00iA /), so 

C3<4|| J D( 2 )y|| 0O ,Mr||/i||+ C3 /2. 

This gives 

sup \\D t Y x+h (t) - D^it)]] = c, < SWD^VW^MrWhl 

0<t<T 

(iv) Set 

Qjf(t) := DiY^\t) - DiY?(t) 
and Q x,h {t) = {Qi'i(t), . . . , Q^n(t))- Using the Taylor expansion we get from ( 13. 14ft . 

D t Y*+ h (s) £ D lm V j {Y*{s))R* rk h {s) ds + 0{\\h\\ 2 ) 

- 1=1 m=l 

+ f D^V 3 iY x (s))-Qf\s)ds 
Jo 

n 

J2D i Y l x+h (s)D^D l V :j (Y x (s)) ■ R x ' h (s) ds + 0(\\h\\ 2 ) 
1=1 

+ f DVVj{Y*{s))-Q°> h {s)ds. 
Jo 

Consider k G {1, . . . , n} and let h = ue^. Define 

< /., (v r • fGfti 

C4 := 04(0;, t, k) := max sup hmsup — nmmt — — 

l<3<iO<Kr V u^O U U 

Note that c 4 is finite because we have \Q x f{t)\ < 8 ||-D^V||oo,m T u for u G (—1,1), by 
(13. 9p . For G (—1, 1) \ {0}, /i = we^ and ft/ = u'e^, (I3.15P implies that, 

9ll^_9llP- = f j2D l Y x+ \s)D^D l V 3 {Y x {s)) ■ ?-^±ds + 0{u) 

ft n TDX,h' 



(3.15) 



ft " T}X,h { a\ 

- \ Y,D i Y l x+h \s)D i ' 1) D l V j {Y x (s)) ■ ; ' ds + 0{u') 
J ° 1=1 ' u ' 

+ f DmVj{Y * {s)) (Qm-9m) ds . 



u u 
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The first two integrals cancel in the limit as -> 0. To see that we can pass to the 
limit, we use the bounded convergence theorem. This theorem is applicable because 
(13. 7ft provides a bound for - R x ' h (s), (13. 8p provides a bound for DiY x+h (s) and we also 
have \\DV>V\\ oo m < oo, by assumption. Letting u,«'-^0 we get 

l imsup 9l£^l _ i iminf 911^1 < T \\ D W V | U M . C4 . 

u^O U u'->0 U 

Since < t < r and j e {1, . . . , n} are arbitrary we see that 

C4<t|| j d( 1 V|| 0O) m-c 4 , 
so C4 = 0; this proves that DikY x (t) exists. The estimate ( 13.101) follows now from ( 13.91) . 
(v) By (I3.15P we get for h = ue k 

; (/> = lin,^" ; : / 



D*Y?(t) = lim 



/ ^ D^D^iY^s)) ■ D k Y x {s) ds 

Jo i=i 

+ / D^V j (Y x (s))D ik Y x (s)ds. 
Jo 



Let x <E K, h <E R n , \\h\\ < 1 and i, k e Put 

C fi:=Cfi(a;,/i,i,A;):= sup || A^+^t) - AfcV a 

0<t<T 

Because of (I3.10p . C5 is finite. For any < t < r and j G {1, . . . , n} we have 

/^y'f/) /; ; ,v/i/) 

E E [D i Yr + \s)D lm V j {Y^ h {s))D k YZ +h {s) 

1=1 m=l 

ft n r -I 

+ / r a^o^^d^w - A^r («))^w d« 
= 1+11. 

We obtain from (j3Tg]l . (I33|) and (fXTUj) . 

i 1 ^ / EE[ I ^rt) ait^w [^^w - w(* 

■'° Z=l m=l 

+ \D l Y l x (s)D k Y x (s) [A ffl ^(s)) - A^'(^(«))] I Ids 
<^ P (2) V1Lm 32 r||/i|| +8||^ 3 )y|U )M 



STATIONARY DISTRIBUTIONS FOR JUMP PROCESSES 25 

as well as 

ft n 

[II] < / £ [ I D l V j {Y*+\s)) [D lk Yr\s) - D tk Y?(s)] \ 
Jo 1=1 L 

+ \D lk Y?(s) [DW+^s)) - D l V j {Y x {s))] 1 1 ds 

<r wdWvw^m-cs + iqwdWvwI 



Combining these two estimates we find for all < t < r and 1 < j < n, 
\D ik Y*+ h {t)- D ik Yf(t)\ 

< 48 WD^VWl M r 2 \\h\\ + 8 WD^VIUmt \\h\\ + r \\D^V\\ oo,M ■ C5. 

Hence, 

c 5 < 48 \\D^V\\l M r 2 \\h\\ + 8 P (3) V||oo,m t \\h\\ + r WD^V^m • c 5 , 
so, recalling (13.61) . 

c 5 < 96 \\D^V\\^ M r 2 \\h\\ + 16 WD^V^m r \\h\\, 

which finishes the proof. □ 

The next step in proving Theorem 13.11 is the following lemma. 

Lemma 3.5. Let g : R n — )■ R be a function in C 2 . Fix < t x < 00 and let Y x {t) 
be the solution of (13.11) . Assume that the explosion time ((x,u) = 00 for all x G R" 
and all u G f2. Let T 4 g 6e defined by (13. 2p . Assume that K C R n , /or ewer?/ t < t\ 
suppTig C FJ and t/iere exists a convex set M C R™ suc/j £/ia£ w) G M /or a// 
x G FJ 3 := Uie-ft: * — ^ an d oj G fi. Assume that ||V||(3),m < 00 anc ^ ^ 



r = A t x := x 



1 

2p( 1 ^||oo,M'^ 1 VO 
Then we have 

(i) For allO<t<f,x£iK and i G {1, . . . , n}, the derivative DiT t g(x) exists and 

(3.16) DiT t g(x) = E A^W) • 

(ii) For allO<t<f,x&K and i, k G {1, . . . t/ie derivative D ik T t g(x) exists 
and 

(3.17) D ik T t g{x) 

= E^D^g(Y x (t)) ■ D lk Y%t) + J2D i Yf(t)D^(D ] g)(Y x (t)) ■ D k Y x {t) 

(iii) For all < t < f and i,k G {1, . . . , n}, i/ie derivative D ik T t g(x) is continuous 
for x E K . 
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Proof, (i) Let < t < f , x e K, fix i e {1, . . . , n} and let ft. = uej. By Taylor's theorem 
and ( 13. 7ft . we get, 



D t T tg (x) = \un Tt9iX + k) ' Tt9{X 



u->0 U 

lim E / g(y* +h (t))-g(y*(t)) 

u->o y it 

lim E f D {1) g{Y*{t)) ■ (y*+fc(g - Y*(t)) 

mo y tt 

, .. ( El<l, m <nDlm9(Z)iY? +h (t) - l?(t))(l£ + *(t) ~ 

+ lim E = — = 

: E (D^g(Y x (t)) ■ DiY*(t)) + lim E (o ( 
: E (D^g(Y x (t)) ■ DiY*(t)) , 



where £ = £ x ,h,t,i,m is an intermediate point between F a: (f) and This yields 



(ii) Fix i,k E {1, . . . , n} and let ft = ue k . We have, using (i), 

n rr , x ,• DiT t g(x + ft) - DiT t g(x) 
D ik T t g(x) = hm 

- lim E /' £>(1) ^( yx+h (*)) • A^ x+h (t) - £> (1) ^(^(t)) • A^ x (t) 

u-tO \ U 

+ Hm E ■ CDWg(y^g)) - DWg(Y*(t))) 
u-¥0 \ u 
= 1 + 11. 



By (13. 9p and bounded convergence theorem, 

I = E (DWg{Y*(t)) ■ D ik Y*{t)) 
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We apply the Taylor theorem, ( 13. 71) and the bounded convergence theorem to see that 



II = lim E 



lim E / E; =1 A^(^ (1) (^g)r W) • (Y x+h (t) - Y'(t))' 



y y u 
= E DiY?(t)D^\D j9 ){Y*(t)) ■ D k Y%t)j . 

This proves (13.171) . 

(iii) By Proposition ^. 31 all derivatives on the right hand side of ( I3.17P are continuous. 
Thus the function D i ] t T t g{x) is continuous for x G K, i, k G {1, . . . , n} and < t < f. 
This proves (iii). □ 



Proof of Theorem VJ.1\ We set 



T '-2\\DWU M ^ 



We will use induction. The induction step is the following. Assume that T s f G C 2 for 
some s G [0, to]- We will show that for all r < r such that s + r < to we have T s+r f e C 2 
and ||T S+J ./||( 2 ) < oo. To show this we use Lemma [331 Put g = T s f and ti = to — s- 
Note that r < r A ti = f and (7 = T s / satisfies the assumptions of Lemma 13.51 Hence 
we obtain that T r+S f = T r g EC 2 . A combination of the estimates (I3.16p . ( I3.17|) . the 
fact that suppT r g C K and the estimates from Proposition 13.31 yield ||T r g||( 2 ) < 00. 

An assumption of Theorem 13.11 states that / G . Hence, T f = f G C 2 . The 
induction step shows that T s f G C 2 for all s < r A t . Subsequent induction steps 
extend this claim to T s f G C 2 for all s < jr A t , j = 2, 3, . . . Therefore, T s f G C 2 for 
all s < t - □ 

Proposition 3.6. Let {T t } t > &e t/ie semigroup given by (12. 6p of the the process (Y t , S t ) 
defined by (12. ip . Let C 2 (IR 2 ) 6e t/ie c/ass of functions given by Definition \2.4\ We have 



T t : L7 2 (E 2 )^C 2 (E 2 ). 

Proof. We will repeat some of the arguments given before the statement of Theorem 
13.11 Note that the SDE (12.11) is of the form (13.11) . By Lemma [2.11 there exists a unique 
solution of (12. ip with explosion time (((y, s), cj) = 00 for all (y, s) G R 2 and u G Q. 
Suppose that / G C 2 (R 2 ). Then supp / C R x [— r, r], for some r > 0. Fix to > 0. By 
(13. 4p . for any t < to, we have 

(3.18) suppTJcK:=Rx [ - r - t || W"^, r + fo||W'||oo] - 
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We have 

K 3 = |J B*((y, 8 ),3) = Rx{-r-t \\W"\\ oo -3,r + t \\W"\\ oo + 3). 

(y,s)GK 

Let 

M = Ex [-r-2t ||W / "||oo-3, r + 2t \\W"\\ oo + 3]. 

By (E3D we have (Y t (y ' s) , S^ s) ) G M for all (y, s) G # 3 . Rewriting (El]) as Q we 
have Vi (2/1,2/2) = (2/1) 2/2, ^2(2/1,2/2) = VK"(yi). Since W E C 5 and since it is periodic, 
we get UVH^Af < 00. Therefore, the solution of (13. 3 j) satisfies the assumptions of 
Theorem 13.11 It follows that for any t <to we have 

T t /GC 2 and ||T t /|| (2) < 00. 

This and (pOID yield T t / G C 2 (E 2 ). □ 

Proof of Proposition\2M Suppose that / G C C 2 (S x E). Then / G C 2 (E 2 ) where / 
is given by Q. By Proposition ESI T t / G C 2 (E 2 ). Using this and OTOj) we get 
T t s / G C C 2 (S x E). □ 
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